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We report the first observation of tlie doubly Cabibbo-suppressed decay, Dt K^ K'^-k~ and 
tiie significantly improved measurement of the doubly Cabibbo-suppressed decay, K^-k^-k^ 
using 605 fb~^ of data collected with the Belle detector at the KEKB asymmetric-energy e+e^ 
collider. The branching ratio with respect to its Cabibbo-favored counterpart are B{Df 
K+K+-K-)/B{Dt K+K--K+) = 0.229 ± 0.028 ± 0.012)% and B{D+ -> K+-k+-k-)/B{D+ 
K~'K'^n^) = (0.569 ± 0.018 ± 0.014)%, where the first uncertainties are statistical and the second is 
systematic. We also report the improved measurement of KgK^ and D'^ Ks''^^ branch- 

ing ratios using the same amount of data samples. The measured branching ratios with respect to 
the Cabibbo-favored modes are B{D+ K%K+)/B{D+ ~+ K%n+) = (0.1899 ± 0.0011 ± 0.0018) 
and B{D+ ^ Kl'K+)/B{Dt K"gK+) = (0.0803 ± 0.0024 ± 0.0012). 



I. INTRODUCTION 

Cabibbo-suppressed (CS) and doubly Cabibbo- 
suppressed (DCS) decays play an important role in 
studies of charmed hadron dynamics. CS decays of 
nearly all the charmed hadrons have been observed, 
while DCS decays have been observed for only 
and D'^ mesons. The naive expectation for the DCS 
decay rate is of the order of tan^ 6c, where Oq is the 
Cabibbo mixing ang le [H, or about 0.29% relative to 
its Cabibbo-favored (CF) counterpart. p7| 

One expects that the branching ratio of 13+ 
K^n^TT~ is about 2ia,n^ 9c relative to its CF coun- 
terpart since the phase space for — > K^tt^tt^ is 
suppressed due to the two identical pious in the fi- 
nal state. [l8| This expectation is consistent with cur- 
rent experimental results Therefore, we also ex- 
pect the branching ratio of D'^ K^K^n~ is about 
l/2tan^6'c relative to its CF counterpart. Lipkin Q 
argues that SU(3) flavor symmetry |19| implies 



B{D+ K+K+TT-) B{D+ K+TT+TT-) 
B{D+ K+K-TT+) B{D+ K-7T+Tr+) 



tan^ Oc 



a destructive interference between color-favored and 
color-suppressed amplitude in ^ K tt"*" Q . How- 
ever, converting experimental measurements includ- 
ing Kg branching ratios to those of involving K is not 
straightforward due to the fact that one must take into 
account the interference between DCS decay and CF 
decay modes where the interference phase is unknown 

0. In Df decays to K and i^°7r+ final states, 
the ratio of CS decay to that of the corresponding CF 

decay may be larger than tan^ 6c , since D'^ K 
decay mode is a CF but color-suppressed decay mode. 
Precise branching ratio measurements of CS and CF 
charm meson decay modes can thus help in under- 
standing the underlying dynamics of these decays. 

The data used in the analysis were recorded at the 
T(4S') resonance with the Belle detector ^ at the 
e^e~ asymmetric-energy collider KEKB [9|. It corre- 
sponds to an integrated luminosity of 605 fb~^. 



II. FIRST OBSERVATION OF DCS DECAY 

IN D+ K+K+TT- AND IMPROVED 
MEASUREMENT OF D+ K+tt+tt- 



where differences in the phase space for CF and DCS 
decay modes cancel in the ratios. The above relation 
does not take into account possible SU(3) breaking 
effects that could arise due to resonant intermediate 
states in the three-body final states considered here 

In addition, decays of charmed mesons play impor- 
tant role in understanding the sources of the SU(3) 
flavor symmetry breaking effects Such break- 

ing effects can originate from strong final-state in- 
teractions or interference between same final states. 
In particular, D+ K°K+ and Df K°tv+ are 
CS decays with the color-favored tree, annihilation 
and penguin diagrams. For decays, the branch- 
ing ratio B{D+ TfK+)/B{D+ if 7r+) devi- 
ates from the naive expectation of tan^ 6c 2] , due to 



We require that the charged tracks originate from 
the vicinity of the interaction point with impact pa- 
rameters in the beam direction (z-axis) and perpen- 
dicular to it of less than 4 cm and 2 cm, respectively. 
All charged tracks are required to have at least two as- 
sociated hits in the silicon vertex detector [l3| , both in 
the z and radial directions, to assure good spatial res- 
olution on the D mesons' decay vertices. The decay 
vertex is formed by fitting the three charged tracks 
to a common vertex and requiring a confidence level 
(C.L.) greater than 0.1%. Charged kaons and pions 
are identified requiring the ratio of particle identifi- 
cation likelihoods, Cx/il^K + ^tt), constructed using 
information from the central drift chamber, time-of- 
flight counters, and aerogel Cherenkov counters [Tl| . 
to be larger or smaller than 0.6, respectively. In ad- 
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dition, we require that the scaled momentum of the 
charmed meson candidate Xp — p* / \/Q.2bE'^^ — 
be greater than 0.5 to suppress combinatorial back- 
ground as well as D mesons produced in B meson de- 
cays. Here p* is the charmed meson momentum and 
i?CM is the total e"'"e~ collision energy calculated in 
the center-of-mass frame. 

In order to suppress the high background level, we 
apply further selection criteria, which are optimized 
using real data samples since there are some discrep- 
ancies between the Monte Carlo (MC) simulation ^12] 
and the data in the relevant distributions. We use 10% 
of the data samples for optimization and the remain- 
ing 90% for the measurement to avoid a possible bias 
when the same samples are used for both optimization 
and the measurement. Assuming no signal in the DCS 
decay channel, we maximize Ns / V-^b, where A/g is 
the CF signal yield which has similar properties to the 
DCS signal and JVb is the background yield from the 
sideband regions in the DCS sample. 

One of the selections related to the finite lifetime of 
charmed hadrons is the reduced (x^/d.o.f) for the 
hypothesis that the candidate tracks for the charmed 
meson decay products arise from the primary vertex. 
The primary vertex is obtained as the most probable 
point of intersection of the meson's momentum vec- 
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FIG. 1: Distributions of M{K'tt+tv+) (top) and 
M(is:+7r+7r") (bottom). The K~Tv+n+ distribution is 
shown on a semilogarithmic scale to make the small back- 
ground visible. Points with error bars show the data and 
histograms show the results of the fits described in the 
text. Signal, Df — + K'^K^tv^ background, and random 
combinatorial background components are also shown. 



tor and the e+e^ interaction region. Because of the 
finite lifetime of Z)+ and mesons their daughter 
tracks are not likely to be compatible with the pri- 
mary vertex. The second requirement uses the angle 
between the charmed meson momentum vector, as re- 
constructed from the daughter tracks, and the vector 
joining its production and decay vertices. In an ideal 
case without resolution the two vectors would be par- 
allel for the signal. The reduced is required to be 
greater than 25 for and 5 for I?+ candidates and 
the angle is required to be less than 1° for and 2° 
for candidates. Tighter requirements on charged 
kaon identification (> 0.9) and Xp (> 0.7) are also 
chosen for the final selection, which improves the sig- 
nal sensitivity. After the additional and tighter selec- 
tion requirements described above, 9.57% of £)+ and 
10.71% of D+ CF signal and 0.06% of D+ and 0.24% 
of DCS background events are retained. In order 
to minimize systematic effects we choose the same se- 
lection criteria for both DCS and CF decay channels. 
The Ktttt and KKtt invariant mass distributions after 
the final selections are shown in Figures [T] and [2] to- 
gether with signal and background parametrizations. 
A clear signal is observed in both DCS decay mass 
distributions. 

The i^TTTT and KKtt invariant mass distributions 
are fitted using the binned maximum likelihood 
method. In all cases the signal probability den- 
sity function (PDF) is parametrized using two Gaus- 
sians with the same central value. Because of K/t: 
misidentification the following reflection backgrounds 
appear in the mass distributions. In — > K~T:'^'n~^ 
(CF) and AT+tt+tt" (DCS) decays there is a contri- 
bution from misidentified K^K~tt^ decays; 
in D'f K^K^TT^ (DCS) decay there is a con- 
tribution from misidentified D'^ K~^Tr~^'!T~ decay. 
The PDFs for the refection backgrounds are deter- 
mined from real data by assigning the nominal pion 
(kaon) mass to a kaon (pion) track. The magnitude 
of each of the reflection background contributions is a 
free parameter in the fit. For the DCS channel, 
the Z?"*" — > Ar"'"7r"''7r~ contribution is not incorporated 
in the fit since it is not significant, but its effect is in- 
cluded as a systematic uncertainty due to fitting listed 
in Table El The D*+ contributions {D*+ D°7r+ 
with D'^ K'^K^) in the CF D+ channel is also in- 
corporated in the CF Df fit as an independent Gaus- 
sian component. A linear function is used for the ran- 
dom combinatorial background for all channels. All 
signal and background parameters for the CF channels 
are fioated. For the DCS channels the mass, width, 
and ratio of the two signal Gaussians are fixed to the 
values obtained from the fits to distributions of CF 
decays. Signal and background yields are left free in 
the fit. From the results of the fits, shown in Figures 
[U and [21 we extract the signal yield for each chan- 
nel, listed together with the corresponding branching 
ratios in Table HI 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



3 




1.94 1.95 1.96 1.97 1.98 1.99 2 2.01 2.02 2.03 

M(K+K"7t+) (GeV/c^) 



> 



c 
o 
> 
m 




■P93 1.94 1.95 1.96 1.97 1.98 1.99 2 2.01 2.02 2.03 

M{K*K*n) (GeV/c^) 



FIG. 2: Distribution of M{K+K-tv+) (top) and 
M{K^ K'^TT~) (bottom). Points with error bars show the 
data and histograms show the resuhs of the fits described 
in the text. Signal, D*'^ background {D'"^ D^tt^ where 
D° K'^K~), D"*" — > K~tt'^tt'^ background, and random 
combinatorial background components are also shown. 



The statistical significance of the Df K^K^tt^ 
signal is calculated using ~'2\n(Cb/ Cg+t) where Ct 
and Cs+b SLie the likelihood values of the fit, without 
and with the signal PDF included, respectively. We 
find — 2 ln(£b/£s+f,) = 83.2 with 1 degree of freedom 
used to describe the DCS signal yield; we obtain a 
statistical significance corresponding to 9.1 standard 
deviations. 

In addition to the backgrounds mentioned above 
there is also the possibility of double misidentification 
leading to contributions from CF events to the DCS 
sample. MC simulation shows that such a contribu- 
tion is flat in the invariant mass distribution and is 
hence included in the combinatorial background de- 
scription. 



The final states in this study have resonant sub- 
structure that can affect the reconstruction efficiency. 
The resonances are relatively well known for the decay 
modes other than Df K^K^tt~ . We used a co- 
herent mixture of resonant contributions according to 
[T^ to generate D+ if ~7r+7r+ decays and calculate 
the reconstruction efficiency. For the D"*" — ^ A'+tt+tt" 
and K^K~it^ decays we used an incoher- 

ent mixture of intermediate states 0. Subsequently 
we varied the contributions of individual intermedi- 
ate states in a correlated manner, within the uncer- 
tainties of the measured branching fractions. The effi- 
ciency calculated from the modified MC sample differs 
from the original one by 1.5% and 2.0% for the 
K~^7r~^TT~ and Df K~^K~Tr'^ decays, respectively, 
and the difference was included in the systematic un- 
certainty of the result. Df K'^K'^Tr~ decays 
were generated according to phase space. For com- 
parison, signal events were generated assuming either 
K*°{S,02)K+ or K*°{1A30)K+ intermediate states. 
The largest relative difference in the efficiency (2.4%) 
was included as a part of the systematic uncertainty. 
Ratios of reconstruction efficiencies for DCS and CF 
decays are found to be 1.042 ± 0.008 ± 0.016 and 
0.963 ± 0.010 ± 0.030 for D+ and D+ decays, respec- 
tively, where the first uncertainty is due to the finite 
MC simulation statistics and the second is the uncer- 
tainty in the resonant structure of the final states. 

With the efficiencies estimated above, we mea- 
sure the inclusive branching ratios of DCS de- 
cays relative to their CF counterparts summa- 
rized in Table [J The product of the branch- 
ing ratios for the two DCS decay modes is found 

0.21) tan^ 9c, where the error is the total uncertainty. 

Several sources of systematic uncertainty cancel in 
the branching ratio calculation due to the similar kine- 
matics of CF and DCS decays (for example, uncertain- 
ties in the tracking efficiencies and particle identifica- 
tion, since the momenta of the final state tracks are 
almost identical). The stability of the branching ra- 
tios against the variation of the selection criteria was 
studied and we observed no changes greater than the 



TABLE I: Measured branching ratios. Brei is the branch- 
ing ratio relative to D'^ K^-r^-k^ for the modes 
and D'^ — > K~-k^ for the D'^ modes. The uncertain- 
ties in the branching ratios are statistical and systematic. 

Decay Mode Msignai Brei (%) 

D+ K+-K+TT- 2637.7 ± 84.4 0.569 ± 0.018 ± 0.014 
D+ K--K+TV+ 482702 ± 727 100 
Dt ^ K+K+TT- 281.4 ±33.8 0.229 ± 0.028 ± 0.012 
D+ A'+Jf"7r+ 118127 ±452 100 



TABLE II: Relative systematic uncertainties in percent, 
where iyBr^i{D+) ^-nd i(d+) systematic uncertain- 
ties for the branching ratio of D"*" and D'^ DCS decays 
relative to their CF counterparts. 



Source 


0'8,,,(D+) 




Fitting 


1.9 


4.2 


MC statistics 


0.8 


1.0 


Reconstruction efficiency 


1.5 


3.1 


Total 


2.5 


5.3 
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TABLE III: Absolute branching fraction for each decay mode and comparisons with previous measurements. The first 
uncertainties shown in the second column are the total uncertainties of our results, and the second are the uncertainties 
in the average CF branching fractions used for normalization [2[. 



Branching fraction 


Belle 


World average [2] 


B[Dt K+K+TT" 


(5.2 ±0.2 ±0.1) 
(1.3 ±0.2 ±0.1) 


X 10^'* (6.2 ±0.7) X 10"" 
X 10"" (2.9 ± 1.1) X 10"" 



expected statistical fluctuations. The systematic un- 
certainties due to the variation of the fit parameters 
are 1.9% for _D+ and 4.2% for D+ branching ratios 
measurements, respectively. Table |TT] summarizes the 
systematic uncertainties in the measurements of the 
branching ratios. 

Using the world average values B{D^ 
K-TT+TT+) = (9.22 ± 0.21)% and B{D+ -> 
K+K^TT+) = (5.50 ± 0.28)% 0, we obtain the ab- 
solute branching fraction for each DCS decay chan- 
nel. Table mil shows the comparison between previous 
results and this work, published in 



TABLE IV; The extracted signal yields (preliminary) from 
the fits to data and corresponding signal efficiencies (e) 
from the simulated events of signal modes. The uncer- 
tainties are statistical only. 

Decay modes Yields e (%) 

D+ ^ K°gK+ 100855 ± 651 12.69 ± 0.01 

D+ K"sK+ 204093 ± 768 13.53 ± 0.01 

D+ K°c;Tr+ 566285 ± 1162 14.19 ± 0.01 

Dt ^ K°gTr+ 17583 ± 481 15.35 ± 0.01 



III. PRECISE MEASUREMENT OF 

D+ K%K+ AND Dt K'-i;-K+ BRANCHING 
RATIOS 

We require similar selection criteria described in 
previous section. In addition, we require that Kg can- 
didates have the invariant mass within ±9 MeV/c^ 
from its nominal Kg mass ^21. Further requirements 
are also imposed to improve the quality of Kg candi- 
dates [HI. and Df candidates are reconstructed 
using a Kg and a charged track. The decay vertex is 
formed by fitting the Kg and the track to a common 
vertex and requiring a confidence level greater than 
0.1%. In order to remove peaking backgrounds from 
— !■ 7r~''7r"*"7r" and K^tt'^tt" decay modes, we com- 
pute the reduced of the vertex assuming that two 
pions from Kg and the charm daughter track arises 
from a single vertex. We require the reduced to be 
greater than 10. To completely remove and 
mesons produced in B meson decays, we require p* to 
be greater than 2.6 GeV/c. Reconstruction efficien- 
cies are found to be 16.6% for _D+ and 18.0% for Df 
in K^K+ final states, and 20.6% for D+ and 22.4% 
for D+ in KgTT^ final states, at this stage. 

Highly asymmetrical Kgh~^ pairs that have in- 
variant mass close to Z)^^ mass region are more 
likely to be background than signal. The selec- 
tion optimization is made on asymmetry, A = 

Pk° Ph+ I Pk° + Ph+ to reject background can- 
didates, where momentum is calculated in the labo- 
ratory frame and refers to either K^ or 7r+. We 
use the off-resonance data sample for optimization and 
the on-resonance data for the measurements as well as 



described in previous section. 

We perform a tuning of MC simulated events [l2| 
intended mainly for the parametrization of the peak- 
ing background under the signal. This background is 
a consequence of particle misidentification and is to be 
discussed in more details later. The tt"*" (K^) momen- 
tum and its resolution are tuned with Z?" —^ tt'^tt^ 
[D^ K^K^) data sample. For K^ momentum 
scale and resolution tuning, — > A'^tt^ data sam- 
ple is used. The tuning method is applied to K^K'^ 
final state simulated data and compared to the real 
data. The agreement between the simulated and the 
real data is significantly improved. 

In the branching ratio measurements, there is a 
peaking background due to particle misidentification. 
In KgK'^ mass window, a peaking struc- 

ture appears from —^ Kgn^ decays when tt"*" 
is misidentified as K^ . A similar peaking structure 
in KgTT^ mass window appears due to the 

misidentification in Df KgK'^ decays. The shapes 
and the yields of these peaking backgrounds are ob- 
tained from the simulated event samples. The sim- 
ulated shape and normalization of the peaking back- 
ground is checked by comparison to the invariant mass 
distributions of selected KgK'^ (Kgir^) pairs with 
K~^{7T~^) mass swapped by the Tr'^{K'^) mass. The 
comparison shows that the simulated peaking back- 
ground of tuned sample correctly describes these com- 
ponents and that the mentioned misidentification is 
indeed the only contribution above the structureless 
combinatorial background. Possible uncertainties in 
the misidentification probabilities are considered as a 
source of systematic uncertainty. 

The KgK~^ and Kgir'^ invariant mass distributions 
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FIG. 3; Invariant mass distribution of selected KgK^ 
(top) and KgTT'^ (bottom) pairs. Points with error bars 
show the data and histograms show the results of the fits 
described in the text. Signals, —* K'^tt'^ (top) and 
D'^ KgK'^ (bottom) backgrounds, and random combi- 
natorial background components are also shown. Inset of 
bottom is an enlarged view of region. 



after the final selections are shown in Figure [3] with 
signal and background parametrizations. Enhanced 
signals of CF and CS decays are observed in both 
distributions. The KgK^ and Kgir^ invariant mass 
distributions are fitted using the binned maximum 
likelihood method. In all cases the signal PDF is 
parametrized using two Gaussians with a common 
mean value. For D+ Kgir'^, we fix the ratio 
of widths and the fractional yields in two Gaussians 
fits due to a lower signal sensitivity. The values of 



the ratio and the fraction are obtained from the fit 
to L)+ ^s^^ mode and are consistent with fit 
results from the simulated events of signal modes. 

/t:^ misidentified backgrounds appear in the way 
described earlier and their normalizations are fixed 
during the fits to mass distributions. Table IIVI sum- 
marizes the extracted signal yields from the fits to 
data and corresponding signal efficiencies from the 
simulated events of signal modes. 

Various contributions to the systematic uncertain- 
ties for the branching ratio measurements are summa- 
rized in Table |Vl The efficiency differences in parti- 
cle identification are estimated and the central values 
of our ratio measurements are corrected. Uncertain- 
ties of such corrections are included in the systemat- 
ics estimate and are found to be 0.90% of the mea- 
sured ratios. We refit data with varying bin size of 
histograms, different fit interval, and changes in the 
combinatorial background PDF. We estimate 0.24% 
and 0.87% of the measured ratio as systematic uncer- 
tainties due to variations in fit methods for and 
modes, respectively. Particle identification and 
its normalization effects of the /n^ misidentified 
background yields in fits are also estimated and found 
to be 0.16% and 0.62% of the measured ratio for 
and Df modes, respectively. Finally, systematic ef- 
fects due to extra constraints in the Df — » Kgir^ 
signal PDF are estimated by re-fitting the data al- 
lowing the fixed parameters to change within their 
one standard deviation. This gives a systematic effect 
corresponding to 0.37% of the measured ratio in Df 
decay modes. 

With the signal efficiencies and the corrections due 
to particle identification efficiency differences, we find 
branching ratios to be B{D+ K'^K+)/B{D+ 
X^7r+) = 0.1899 ± 0.0011 ± 0.0018 and B{D+ 
K+TT+)/B{Df K+K+) = 0.0803 ±0.0024 ±0.0012 
where the first uncertainties are statistical and the 
second are systematic. These are the most precise 
measurements up to date and are compared to the 
present world average values in Table IVII 



TABLE V: Relative systematic uncertainties in percent 
(preliminary), where o-^(£)+) and cr^j^o+j are systematic 
uncertainties for branching ratios of D"*" and Df decays. 
Sources include particle identification (PID), fit methods, 
peaking background, and Df signal PDF. 



Source 


^R(D+) (%) 




PID 


0.90 


0.90 


Fit methods 


0.24 


0.87 


peaking background 


0.16 


0.62 


Df signal PDF 




0.37 


Total 


0.95 


1.45 
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TABLE VI: Branching ratios for _D+ and Dt (prelimi- 
nary), and comparisons with previous measurements. The 
uncertainties shown include the statistical and systematic 
uncertainties of our results. 

Branching Ratio Belle exp. World average [2| 



B{D+^K|7r+) 
B(D+^-K'°K+) 



(19.0 ±0.2)% 
(8.0 ± 0.3)% 



(20.6 ± 1.4)% 
(8.4 ±0.9)% 



IV. CONCLUSIONS 

Using 605fb^^ of data collected with the Belle de- 
tector at the KEKB asymmetric-energy e+e" col- 
lider we have observed for the first time the de- 
cay K^K^n' with a statistical signifi- 
cance of 9.1 standard deviations. This is the first 
DCS decay mode of the meson. The branch- 
ing ratio with respect to the CF decay is (0.229 ± 
0.028 ± 0.012)%, where the first uncertainty is sta- 
tistical and the second is systematic. We have 
also determined the DCS decay branching ra- 
tio, B{D+ K+Tr+TT-)/B{D+ if-7r+7r+) 
(0.569 ± 0.018 ± 0.014)%. with significantly improved 
precision compared to the current world average [![. 
We find the product of the two relative branching ra- 



tios, 



B(DJ^K+K + 7T-) B(£)+^g+7r+7r~) 



B(Dt^K+K-^+) B(D+-,K-^+^ + ) ' ^ be (1.57 ± 

0.21)tan^6'c- This is consistent with SU(3) flavor 
symmetry within 3 standard deviations; note that the 
effect of (different) resonant intermediate states is not 
taken into account in the prediction Q . An amplitude 



analysis on a larger data sample may allow a more pre- 
cise test of SU(3) flavor symmetry to be performed. 

We also have measured branching ratios with 
respect to the CF mode and values are B{D^ 
K°K+)/B{D+ K°n+) = 0.1899 ± 0.0011 ± 0.0018 
and BiD+ K%tt+)/B{D+ K%K+) = 

0.0803 ±0.0024 ±0.0012. Using the world average val- 
ues of CF decay modes we obtain branching frac- 
tion values olB{D+ K%K+) = (2.75±0.08)x 10-3 
and B{D+ -> KItt+) ^ (1.20 ± 0.09) x lO'^. These 
are consistent with the present world average 
and best measurements up to date. The ratio 
B{D+ -> KlK+)/B{D+ Kl'K+) = 2.29 ± 0.18, 
which implies possible flavor SU(3) breaking or 
different flnal-state interaction between and 
decays. 
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